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Abstract 

Superconducting structures with the size of the order of the superconducting 
coherence length £(T) have a critical temperature T c , oscillating as a function 
of the applied perpendicular magnetic field H (or flux <I>). For a thin-wire su- 
perconducting loop, the oscillations in T c are perfectly periodic with H (this 
is the well-known Little-Parks effect), while for a singly connected super- 
conducting disk the oscillations are pseudoperiodic, i.e. the magnetic period 
decreases as H grows. In the present paper, we study the intermediate case; a 
loop made of thick wires. By increasing the size of the opening in the middle, 
the disk- like behaviour of T C (H) with a quasi-linear background (character- 
istic for '3-dimensionaF (3D) behaviour) is shown to evolve into a parabolic 
T C {H) background ('2D'), superimposed with perfectly periodic oscillations. 
Lhe calculations are performed using the linearized Ginzburg-Landau theory, 
with the proper normal/vacuum boundary conditions at both the internal and 
the external interface. Above a certain crossover magnetic flux <£, T C (<I>) of 
the loops becomes quasi-linear, and the flux period matches with the case of 
the filled disk. This dimensional transition is similar to the 2D-3D transition 
for thin films in a parallel magnetic field, where vortices enter the material 
as soon as the film thickness t > 1.8 £(T). For the loops studied here, the 
crossover point appears for w ~ 1.8 £(T) as well, with w the width of the 
wires forming the loop. In the 3D regime, a "giant vortex state" establishes, 
where superconductivity is concentrated near the sample's outer interface. 
The vortex is then localized inside the loop's opening. 

74.60.Ec, 74.25.Dw, 73.23.-b, 74.20.De, 74.76.-w 



Typeset using REVTgX 



1 



I. INTRODUCTION 



The nucleation of superconductivity in mesoscopic samples has received a renewed in- 
terest after the development of nanofabrication techniques, like electron beam lithography. 
A superconductor is called mesoscopic when the sample size is comparable to the magnetic 
penetration depth A(T) and the superconducting coherence length £(T). In the framework 
of the Ginzburg-Landau (GL) theory, £(T) sets the length scale for spatial variations of 
the modulus of the superconducting order parameter |^|. The pioneering, work on meso- 
scopic superconductors was carried out already in 1962 by Little and Parkso, who measured 
the shift of the critical temperature T C (H) of a (multiply connected) thin- walled Sn mi- 
crocylinder (a thin-wire "loop") in an axial magnetic field H. The T C (H) phase boundary 
showed a periodic component, with the magnetic period corresponding to the penetration 
of a superconducting flux quantum $o — h/2e. The Little-Parks oscillations in T C (H) are a 
straightforward consequence of the fluxoid quantization constraint, which was predicted by 
F. London!. This condition can be easily understood by integrating the second GL equation 

J= |uf (v6-^A) = \*\ 2 v s (1) 



V % 

along a closed contour. J = j/(2eh/m*) is here the normalized current density, v s is the 
(normalized) superfluid velocity, and 5 is the phase of the complex order parameter ^ = 
\^\e t5 . This yields: 

J $ r - 

dl + 2ir—=fV5-dl = L2ir £, = ■■■ , —2, —1, 0, 1, 2, • • • (2) 



|*| 2 $o J 

where we used Stokes theorem § A • dl = $, with $ the magnetic flux threading the area 
inside the contour. In other words, when a non-integer magnetic flux $/$o is applied, 
a supercurrent J has to be generated in order to fulfill Eq. (0). The integer number L 
is the phase winding or fluxoid quantum number and gives the number L of flux quanta 
$o penetrating the sample. Since, for a cylindrical geometry, the different L states are 
eigenf unctions of the angular (or orbital) momentum operator as well!, L is often called the 
angular momentum quantum number. The fluxoid quantization constraint (Eq. @) is the 
equivalent of the Bohr-Sommerfeld quantum condition for Cooper pairs!. 

A few years later, Saint- James calculated the T C (H) (or the nucleation magnetic field) of 
a singly-connected cylinder! (a mesoscopic "disk"). Taking into account the analogy with 
the situation of a semi-infinite superconducting slab in contact with vacuum!, the critical 
field was called H c3 (T) in this case, since superconductivity nucleates initially near the 
sample interface. As we will show further in this paper, a substantial enhancement of the 
nucleation field, above the H c3 (T) for a semi-infinite slab, can be obtained in mesoscopic 
samples, where the surface-to-volume ratio is large. In such case, surface effects play the 
role of the bulk effects. Therefore, we will rather use the notation H* 3 (T). For a disk, the 
limiting value H* 3 (T) = H c3 (T) = 1.69 /f c2 (T) is found as $ — > oo (or the radius of the disk 
— > oo). The field H c2 is the upper critical field of a bulk Type II superconductor, at which 
the Abrikosov vortex lattice is formed. 

The T C (H) phase boundary (or H* 3 (T)) of the disk shows, just like for the usual Little- 
Parks effect in a multiply connected sample (loop), an oscillatory behaviour. In this case, 
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the oscillation period of T C (H) is not constant, but slightly decreases as H increases. When 
moving along T C (H), superconductivity concentrates more and more near the sample inter- 
face as H grows. A giant vortex state is formed: a "normal" core carries L flux quanta, and 
the 'effective' loop radius increases, resulting in a decrease of the magnetic oscillation period. 
Here as well, fluxoid quantization (Eq. (0)) is responsible for the oscillations of T c versus H. 
An experimental verification of these predictions was carried out later on by Buisson et a/.i 
and by Moshchalkov et a/.i. 

Currently, manv different sample topologies are studied: superconducting networks!, 
antidot systemst2nl3, samples consisting of sharp corner SB, etc. With the use of sub- 
micron Hall probe microscopylll for measuring the magnetic response of a superconductor, 
it has become possible to probe samples deep in the superconducting stately (i.e. at tem- 
peratures far below T c ). A lot of recent theoretical activities have been devoted to the 
magnetization and vortex configurations in superconducting disks of different size MM. 

Loop structures have also been studied extensively in the past years. A large portion 
of the theoretical research has been focused on the transition between two different fluxoid 
states L — > L + lEiE, but also experiments, among which susceptibility measurements close 
to Tjli, studies of energy gap spectrum using a single-electron transistoiS, and scanning 
Hall probe measurements on ensembles of mesoscopic Al loopsS, were carried out. 

As it has already been mentioned, the two limiting cases, "thin- wire loop" and "disk", 
are well understood as far as the nucleation field H* 3 is concerned, but the intermediate case 
is not. In the early paper from Saint- James and de Gennesi, H* 3 (T) has been calculated 
also for a film exposed to a parallel magnetic field, where surface superconductivity can 
grow along two superconductor /vacuum interfaces. For low magnetic fields, the two surface 
superconducting sheaths overlap, and, as a result, T c versus H becomes parabolic, which 
is characteristic for 2D behaviour. When increasing the field, a crossover to a linear T C (H) 
dependence (3D) occurs at t « 2£(T), with t the film thickness. Shortly after, it was shown 
that vortices start to nucleate in the film at this dimensional crossover point (t = 1.8 £(T))il. 

The goal of the present paper is to study the phase boundary T C (H) of loops made 
of finite width wires (or disks with an opening in the middle, increasing in size). In a 
Type-II material, superconductivity is expected to be enhanced, with respect to the bulk 
(H* 3 (T) > if c2 (T)), both at the external and the internal sample surfaces. As for a film in 
a parallel field, a 2D-3D dimensional crossover can be anticipated, since the loops may be 
simply considered as a film, which is bent such that its ends are joined together. We calculate 
the phase boundary T C (H) as the ground state solution of the linearized first GL equation 
with two super conductor /vacuum interfaces. This calculation has been suggested already 
by several authors0H, but to the best of our knowledge, it has not been carried out so far. 
We will show that the T c ($) of the loops, for low applied magnetic flux (corresponding to 
the 2D regime), can be described within a simple London picture, where the modulus of the 
order parameter is spatially constant. As the flux increases, the background depression 
of T c ($) is changed from parabolic (2D) to quasi-linear (3D), which indicates the formation 
of a giant vortex state, where only a surface sheath close to the sample's outer interface is 
in the superconducting state. Moreover, the oscillation period of T c ($) becomes identical 
for the loops as for the full disk, as soon as the transition to 3D behaviour has taken place. 
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II. THE LINEARIZED GL EQUATION 



The linearized first GL equation to be solved in order to find T C (H) is: 

— (-i/iV - 2eA) 2 m = I - orl * . (3) 
2m* v w 

This equation is formally identical to the Schrodinger equation for a particle with a charge 
2e in a magnetic field. Here, at the onset of superconductivity T « T C (H), the nonlinear 
term /3|\l/| 2 \l/ <C | — a| \I/ can been omitted, and demagnetization effects along the field 
direction do not need to be considered. In this regime, the ^-dependence disappears from 
the equations and therefore an infinitely long cylinder and a disk have identical T C (H) 
boundaries. It is further assumed that the penetration depth A(T) sample size, so that 
HqH = rot A. It is important to note that H is the applied magnetic field, so fields induced 
by supercurrents J (Eq. (Q)) are not taken into account here. The more Type II (the higher 
k) the superconducting material is, the larger range of validity our calculations have. The 
eigenenergies | — a\ can be written as: 

I - I = ^ = (-, _ T_\ U) 

1 2m*£ 2 (T) 2m*£ 2 (0) \ T c0 ) ' U 

T c0 being the critical temperature in zero magnetic field. From the energy eigenvalues of 
Eq. (H), the lowest Landau level | — Qlll(H)\ is directly related to the highest possible 
temperature T C (H), for which superconductivity can exist. 

By varying the topology of the sample ("nanostructuring"), the lowest Landau level 
I — can tuned by confinement of the superconducting condensate. Several examples 

of this concept can be found in Refs.i'11. Indeed, the solution \I/ of Eq. (|3]) has to fulfill the 
Neumann boundary condition 

(-thf - 2eA)v\ = , (5) 

at the sample interfaces b. This requirement guarantees that the supercurrent does not have 
a component perpendicular to a superconductor /vacuum interface. 

For the loop geometries, we choose the cylindrical coordinate system (r, ip) and the 
gauge A = (fi Hr/2) e^, where e v is the tangential unit vector. The exact solution of the 
Hamiltonian (Eq. (|3|)) in cylindrical coordinates takes the following formQ'llilfll 

*($, <p) = e"^ (p\ m exp K(-n, L + 1, $/$ ) (6) 



k $o/ V 2$ 

K(a, c, y) = ci M(a, c, y) + c 2 U (a, c, y) . 

Here $ = fioHrcr 2 is the applied magnetic flux through a circle of radius r. The number 
n determines the energy eigenvalue. Most generally, the function K(a, c, y) can be any 
linear combination of the two confluent hypergeometric functions (or Kummer functions) 
M(a,c,y) and U(a, c, but the sample topology puts a constraint on ci, c 2 , and n, via 
the boundary condition (Eq. (|J)). 

The eigenenergies of Eq. @ (the Landau levels) areB 
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a = (2n + 1) = Tiu [n + - 

1 2m* K J V 2 



(7) 



where to = 2efioH/m* is the cyclotron frequency. The parameter n depends on L and is not 
necessarily an integer number, as we shall see later. With Eq. (f|) this can be rewritten as: 



i 



T C (H)' 



T 
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cO 



1\ $ $ 



(8) 



where $ = fioHnrl is arbitrarily defined. 

The bulk Landau levels can be found when substituting n — 0, 1, 2, • • • in Eqs. (0) and (||), 
meaning that the lowest level n = corresponds to the upper critical field hqH C 2(T) = 
$o/ (27r£ 2 (T)). Let us note that the lowest Landau level (n = 0) for a bulk superconductor is 
degenerate in the phase winding number L, and therefore the eigenfunction can be expanded 
as \& = Interference patterns between the different functions ^ l give rise to 

aperiodic vortex statesH. 

The boundary condition (Eq. (|5])), in cylindrical coordinates, can be simply written as: 



9|*(r) 



Or 



with a superconductor/vacuum interface at a radius r Q . Using ^MlMM 



and dU ^' yS) = —a U(a + l,c+l,y) for the derivatives of the first and second type of Kummer 
functions, respectivelyil, and inserting Eq. 



dy 



(9) 



fM(a + l,c + l,y) 



Cl 
C-2 



D into Eq. gives: 

2n $ 

— M(-n + l,L + 2,$/$ ) 



L - ^ W (-n, L + l, $/$ ) + 2 n -J- ?7 (-n + 1, L + 2, $/$ ) 



+ 
0. 



(10) 



which has to be solved numerically for each integer value of L, resulting in a set of values 
n(L, $), with $ = /i if7rr;L 

For a disk geometryii'll, we have to take C2 = in Eqs. (|6]) and (|T0|) in order to avoid the 
divergency of U(a, c, y — > 0) = oo at the origin. Selecting the lowest Landau level at each 
value $, one ends up with a cusp-like T C (H) phase boundary!, which is composed of values 
n < in Eq. ([/]), thus leading to H* 3 (T) > H c2 (T). A similar calculation was performed for 
a single circular microhole in a plane film ("antidot" )§, where C\ = in Eqs. @ and ([10"D , 
since M(a, c,y —>■ oo) = oo. Here as well, the lowest Landau level consists of solutions with 
n < 0. At each cusp in T c ($), the system makes a transition L — > L± 1, i.e. a vortex enters 
or is removed from the sample. 

The loops we are currently studying have two superconducting/vacuum interfaces, one 
at the outer radius r a , and one at the inner radius ?v Consequently, the boundary condition 
(Eq. (|I0"D) has to be fulfilled at both r G , and ?v As a result, we have a system of two 
equations and two variables n and c 2 (ci = 1 is chosen), which we solved for different values 
of x = ri/r Q . In the rest of the paper we will define $ = fioHTcr^, where r Q always means the 
outer loop radius. 
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III. THE LONDON LIMIT 



The usual description of the Little-Parks effectii is given in terms of the London limit, 
where |^| is spatially constant. This approximation, of course, is valid when the wire, 
forming the loop, is very thin (x ~ 1), or when we define the loop aspect ratio z as 



r„ - r, ; 1 - x 



(11) 



r Q + Ti 1 + x 

this condition implies z C 1. 

The solution of the linearized GL equation (Eq. (^)) become; 

era = e<o) = ( 1 + z )- 2i w + 2i'H— ) • (12) 

with $ m = [loHnr^, r m is the mean radius of the loop. Note that this definition of flux 
$ m is different from <3> in the previous section. The lowest eigenvalues are obtained when L 
is the integer number closest to —2 ($ m /$o) z/hix. We will further compare this equation 
with our more exact results, calculated from the scheme presented in the previous section. 

Since in the original paper by Groff and Parks0, thin-wire loops (z <C 1) were investi- 
gated, the logarithm was expanded in a Taylor series, which gives, up to the order z 2 : 

rl _ r 2 m T C (H)\ _ f L _<^\ 2 + 4 z2 /<KA 2 ^ ^ 



The first term on the right hand side of Eq. ( p!3[ ) is the periodic part of the T c reduction 
(i.e. the Little- Parks effect), while the second term is a monotonic parabolic background, 
which is identical to the Tc{H) expression for a plane film of thickness t = 2 z r m in a parallel 
magnetic field. In Ref.E3, a substitution is performed, which splits the right hand side of 
Eq. ( |l~2|) in a periodic term and a monotonic background T c . The latter becomes: 



£ 2 (T C ) v ' \\ 2z ) \l-z 



which is parabolic with $ m . This substitution is only valid for thin-wire loops, and Eq. Q14D 
transforms of course into the last term of Eq. (|1^) for z< 1. 



IV. RESULTS 



Fig. [T] shows the Landau level scheme (dashed lines) calculated from Eqs. (f§), and fliPf ) 



for loops with a different inner radius x = ri/r a (a) x = 0.1, (b) x = 0.3, (c) x = 0.5, 
(d) x = 0.7, and (e) x = 0.9. The applied magnetic flux $ = fi Hnr 2 is defined with respect 
to the outer sample area. The T C (H) boundary (or | — ozlll(H)\) is composed of $ solutions 
with a different phase winding number L and is drawn as a solid cusp-like line in Fig. [I]. At 
$ ~ 0, the state with L = is formed at T C (H) and one by one, consecutive flux quanta L 
enter the loop as the magnetic field increases. Each state L approximately has a parabolic 
dependence | — a(H)\, close to T C (H). Like for the disk@!, where $ w <3> (L + L 1 / 2 ), we 
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have here as well L <3>o < indicating the overall diamagnetic response of the sample. As x 
increases, the oscillations in T c ($) change from cusp-like, to very pronounced local extrema 
for x = 0.9. In the limit of vanishing wire width (x — > 1), L is the integer closest to Q m /§o, 
and therefore the response of the loop is alternating between diamagnetic and paramagnetic, 
as the flux varies. 

The solid and dotted straight lines in Fig. [I] are the bulk upper critical field H c2 (T) 
and the surface critical field H c3 (T) for a semi-infinite slab, respectively. In these units 
the slopes of the curves (see Eq. ([8])) are e = 2 for H c2 (substitute n = in Eq. (||)) 
and e = 2/1.69 for H c3 . The ratio r\ = e(H c2 ) / e(H c3 ) = 1.69 corresponds then to the 
enhancement factor H c3 (T) / H c2 (T) at a constant temperature. For the loops we are studying 
here r\ = e(H c2 ) / e(H* 3 ) is no longer a constant, but varies with the magnetic field. 

The energy levels below the H c2 line (solid straight line in Fig. [I]) could be found by 
fixing a certain L, and solving Eq. (|10"D for a small $, until a set (n, c 2 ) is found with n < 0. 
These values were always put in as starting values for a slightly higher $. A trivial solution 
of Eq. ( f[0|) is obtained for n = 0, 1, 2, • • • . Both confluent hypergeometric functions reduce 
to M(-N,L + l,®/$ ) = 1 and U(—N, L + l, $/$o) = 1, and thus c 2 = -c x . Inserting this 
into Eq. flfj) gives \l/($, (p) — everywhere. With this method, we were able to find solutions 
with n < numerically. Note that the lowest Landau level always has a lower energy 
| — a($)| than for a semi-infinite superconducting slab, which implies H* 3 (T) > H c3 (T). 

The dash-dotted curve in Fig. [I] gives the result obtained with the London limited (see 
Eq. (|L2"D). In Fig. |l|a {x = 0.1) the deviation from the exact solution of Eq. (|) appears 
already for L > 1. For low flux, the result from Eq. ([12]), has clearly higher energy \ —az\ than 
the surface critical field H c3 (T). At $ ~ 7$o, this curve even crosses the bulk H c2 (T) line, 
which is clearly unphysical. For x = 0.3 (Fig. [l|b) the London limit is valid up to $ ~ 4$ ? 
up to $ ~ 8 $o f° r x — 0-5 (Fig. |I|c), and for x = 0.7 (Fig. |I]d) it is a good approximation in 
the whole flux interval of our calculations. Finally, for x = 0.9 (Fig. |l|e) the assumption of a 
spatially constant \^\ gives a T c ($), which can not be distinguished from the exact solution 
of Eq. (|). 

In Fig. 0a (x = 0.1) the background depression of T c is quasi-linear, just like for the case of 
a filled disk, for x = 0.3 (Fig. |I]b) T C (H) has a rather parabolic background for flux $ < 7 $ 0) 
and becomes quasi-linear at higher flux. This becomes more clear for x = 0.5 (Fig. [l|c). 
Here, the crossover point from parabolic to quasi-linear appears at about $ m 14 $o- In 
Figs. |l|d and |l|e finally (x = 0.7, x = 0.9), the background is parabolic in the entire flux 
regime and can be very accurately described by Eq. (|T^). Simultaneously, as x increases, 
the cusps in T C (H) become more and more pronounced, until the usual Little-Parks effect is 
recovered for x = 0.9 (see Eq. (|T3|)), where sharp local minima and clear maxima in T C (H) 
are seen. 

V. DISCUSSION 

In order to study the spatial variation of the order parameter |^|, we will use Abrikosov's 
definition of the flatness parameter (3a =< |^| 4 > / < l^) 2 > 2 , where the brackets indicate 
the average over the actual sample are (not including the middle opening). Abrikosov in- 
troduced this parameter (3a in order to find the most favorable order parameter distribution, 
in a bulk system, near the H c2 line. In this language, (3a = 1 means a completely flat profile 



7 



of \^f\, corresponding to the London limit. For comparison, we mention that (3a = 1.16 for 
a triangular Abrikosov vortex lattice. 

In Fig. we plot the modulus of the order parameter \^> \ for the case x = 0.1, at $ = 9 $oi 
for the states a) L = 2, b) L = 3, c) L = 4, d) L = 5, and e) L = 6. The |\&| values have 
been normalized to 1 at the sample's outer interface r = r a . The dark area is the region 
outside the sample. The maximum in |^| for L = 2 lies at r ~ (r* + r G )/2, and gradually 
shifts to the outer sample edge as L — > 6, which is the ground state solution of Eq. (|3|). The 
spatial modulation of |^| is considerable ((3a > 1.16) for all \^\ patterns shown here. At 
T = T C (<I>) (Fig. (L = 6, (3a — 1.59) the sample is in the "giant vortex state", with a 
normal core containing 6 flux quanta $0; an d a surface superconducting sheath at the outer 
sample edge. 

The |\&| profiles for a loop with a larger inner radius (x = 0.5) are shown in Fig. |3|. Here as 
well, we have chosen $ = 9 $o> an d the same normalization |\I>(r )| = 1. a) L = 2, b) L = 3, 
c) L = 4, d) L = 5, and e) L = 6. For L = 2 (Fig. |3]a), |^| has a maximum at r = ?y 
For higher L the order parameter distribution flattens until it reaches (3a = 1-06 for L = 4. 
Then, for the ground state energy (L = 5) (Fig. |3|d), the maximum in | \I/ 1 moves outward, 
but we should notice that for this state (3a = 1-01 only, which means that superconductivity 
nucleates in a quasi-uniform way. Indeed, the exact T c ($) and the London limit result are 
still very close to each other at $ = 9 $o ( see Fig- ^p)- Although multiple flux quanta L 
are threading the middle opening of the loop, we can not, in the strict sense, speak about a 
giant vortex state here. First of all, there is no real 'normal core' within the sample area, 
and secondly we are not dealing with a surface superconducting state in this case. The 
surface-to-volume ratio is so large that the whole sample area becomes superconducting at 
once. On the contrary, for the disk, strong spatial gradients of are responsible for the 
spontaneous breaking of superconductivity in the giant vortex core, while only a surface 
sheath is superconducting. 

It is worth noting that for all of the states shown in Fig. |2| (x = 0.1), the width of the 
wire w = r a — ri = (1 — x) r Q > 2.7£(T) (see also Fig. |l]a), while for the loop with x = 0.5 
the different L states from Fig. |3] correspond to w < 2.1 £(T) (see also Fig. [l|c). At this 
point we want to remind that in a thin film of thickness t in a parallel field H, a dimensional 
crossover is found at t = 1.84 £(T). For low fields (high £) T C (H) is parabolic (2D), and 
for higher fields vortices start penetrating the film and consequently T C (H) becomes linear 
(3D)E3. In Figs. |l]a, b, and c the small arrow indicates the point on the phase diagram 
T c ($) where w = 1.84£(T). For loops with larger x this point lies outside the flux regime of 
the calculations. For the loops as well, the dimensional transition shows up approximately 
at this point, although the vortices are not penetrating the sample area in the 3D regime. 
Instead, the middle loop opening contains an integer number of flux quanta L$o- in the 
present case, the 2D-3D crossover roughly occurs at the value $ where the London limit 
result (Eq. (|T2"|)) fails to describe the exact T c ($). 

In order to compare the flux periodicity of T c ($), we have replotted, in Fig. |], the lowest 
energy levels of Fig. |I] as rj" 1 = e(H* 3 )/e(H C 2), which is the inverse enhancement factor at 
a constant temperature. In this representation, the dotted horizontal line at rj^ 1 = 0.59 
corresponds to the surface critical field line H c3 (T). The nucleation field of a disk H* 3 (T) > 
1.69 H C 2(T) (i.e. 77 > 1.69), and for a circular microhole in an infinite film ("antidot" )§ 
#* 3 (T) < 1.69 H c2 (T) (77 < 1.69). As $ grows (the radius goes to infinity) the H* 3 (T) of 
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both the disk and the antidot approaches the H c3 (T) line. 

At this point, we would like to come back to the explanation why the notation H* 3 (T) 
was used for the nucleation field. For a thin film (thickness t) with the field H applied 
parallel to the surfaces, T C (H) is found from the second term of Eq. ([13]). With t = 2zr m 
and fioH C 2 = ^o/(2 7r^ 2 (T)), this becomes 

KziT) = ^^ T) H c2 {T) = 2.04 ^p- H c3 (T) , (15) 

which implies directly H* 3 (T) > H c3 (T) for very thin films t < 2.04f(T)H. Note that 
t = w ~ 2.04 £(T) is very close to the 2D-3D crossover point in films and in loops. Of 
course, the possibility for nucleation fields H* 3 > H c3 is not very special, but it still creates 
a lot of confusion. In thin films, for example, the critical field H* 3 (T) is often denoted 
as H C 2 t j_ in a perpendicular magnetic field, and H&n in a parallel field, which would mean 
# c2j ||(T) > H c3 {T) for t < 2.04 f (T) (Eq. (0)). Van Gelder studied the H* 3 of a semi-infinite 
film which is bent over a certain angle (a "wedge") and even called it H ci , since it exceeds 
H c3 at small angles!!. Since this is just a finite size effect and indeed no new nucleation 
mechanism is involved, the existence of a H c/ ± > H c3 was called a misinterpretation by FinkH. 
For all these reasons, it is safe to use the notation H* 3 for the nucleation field. Enhancing 
the nucleation field H* 3 {T) can be realized by confining the superconducting condensate in 
a mesoscopic sample. The smaller the sample size, the lower the Landau level | — ollll{H)\ 
becomes, when Neumann boundary conditions (Eq. (|^)) are imposed. The idea of having 
H* 3 > H c3 = 1.69 H c2 was used to explain anomalously high values of H* 3 and attribute this 
to the existence of small surface irregularities (of the order of £), impurities, grain boundaries 
or concentration gradients near the interfaced. 

For all the loops we study here, the presence of the outer sample interface automatically 
implies that H* 3 (T) > H c3 (T) is enhanced (rj > 1.69), with respect to the case of a flat 
superconductor-vacuum interface. For loops with a small x, the T c ($) boundary very rapidly 
collapses with the T c ($) of the dot (for $ > 4$o) (v becomes the same). Since both the 
flux periodicity and the background depression of T c become identical to the disk, a giant 
vortex state can be anticipated for the loop as well. The presence of the opening in the 
sample is not relevant for the giant vortex formation in the high flux regime. Indeed, for 
the ground state level of Fig. [2] (x = 0.1, L = 6), which is marked with a circle in Fig. |], 
the giant vortex state has clearly been formed already. On the contrary, for x = 0.5, L — 5, 
at $ = 9 $o, marked with a square in Fig. [|, superconductivity nucleates in a uniform way 
(Pa = 1.01) (see Fig. 0). 

The formation of the giant vortex state at high $ can be understood when writing the 
GL free energyES 



h 2 

F = F n + S < a\^\ 2 >+S — < 

2m* 



-zV - —A ) 
$ 



>+•■■, (16) 



with the sum of the last two terms vanishing at the phase boundary T C (H), where ^ — > 0. 
F n is the free energy of the system in the normal state. The notation <> denotes the average 
over the sample area S, and \a\ = h 2 /2m*£ 2 (T). This yields (see Eq. ([[])): 



em 



< 




2 v 2 > + < (V|tf | 


) 2 > 


< 




2 > 



(17) 
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The solutions from Eq. (|6|) which fulfill the boundary conditions (Eq. @) in r = ri and 
r = r Q need to be inserted in this equation, and v s is determined from Eq. (y). A fast 
calculation shows that the relative contribution from spatial gradients (the second term in 
the numerator of Eq. (fL7|)) is 47 % for the lowest level of Fig. ||e (x = 0.1, L = 6), while it 
is less than 4 % for the ground state of Fig. [|d (x = 0.5, L — 5). The energetically most 
favourable balance between these two contributions is strongly affected by the boundary 
conditions (Eq. |5|). In thin-wire loops, for example, bending of \^/\ on a scale smaller than 

the coherence length £ will result in a large contribution of < ^V|\I/|) >. Therefore, the T c 
reduction is only determined by the averaged supercurrent kinetic energy < |\I/| 2 i^ 2 > in 
this case, which is equivalent to the London limit. 

An examination of the supercurrent flow profiles shows that, as in the case of a filled 
disk, there is a paramagnetic contribution close to the sample center for L ^ 0, while near 
the outer sample interface diamagnetic currents are flowing. The supercurrent density J can 
be obtained by inserting the general solution (Eq. @) in the second GL equation (Eq. (|[)): 

J(L,^) = \^\ 2 v-^^L-^j^ L 2 exp(-|-) [K(-N,L + l,^/%)] 2 e v (18) 

At a radius r, corresponding to integer flux quanta $ = L <3> the current orientation changes 
sign. It can happen, however, that the switching radius r < r«, lies inside the middle opening. 
For the state L = there are only diamagnetic currents. In contrast to this, for example, for 
a loop with x = 0.5, the states L = 2 (Fig. |3|a) and L = 3 (Fig. |3]b) do not carry paramagnetic 
currents at $ = 9 $o- At the ground state level (L = 5) (Fig. |3]d), which is very close to the 
London limit solution, J oc v s oc (L — $/$ )/ r (Eq. (0)). This is shown in the right inset 
of Fig. f|, where the supercurrent has been normalized to —1 at the outer radius r Q . For 
the loop with x — 0.1 (see Fig. |2]e), in the ground state L = 6, the dominating diamagnetic 
supercurrents are flowing in the vicinity of r = r D , the currents are paramagnetic for slightly 
lower r, and J is vanishing in the core of the sample < r < r G /3 (left inset of Fig. [|). 

As a last point, we discuss now the periodicity A<3> of T c ($), for the different L states 
in the loops. These results are shown in Fig. g. For the disk, the first cusp in T c ($), where 
L goes from — ► 1, occurs at $ = 1.92 $ , so A$(L = 0) = 3.85 $ . The period A<3> goes 
down for increasing L, until it reaches the asymptotic limitl: 

A$ = $ (l + (2 v $/$or 1/2 ) (19) 

Since r] is decreasing with $, A$ is a weakly decreasing function of L, at high <3>. The symbols 
in Fig. |^ correspond to the period A$ in units of $o for the different L states. Since at 
the highest $, T c ($) is calculated in steps of ~ 0.07 $o, the error on A$ is of this order. 
The filled squares are the periods for the filled disk. The periodicity of T c ($) in the case of 
loops behaves differently: for x = 0.1 (filled up-triangles in Fig. |5|) A$(L = 1) = 1.82 $ is 
larger than for the filled disk, then A$(L = 2) jumps below the corresponding value for the 
disk, and for higher L > 3, the period A$ matches with the disk behaviour. Consequently, 
the giant vortex state builds up for L > 3, i.e. $ ^ 5 $o- A similar analysis can be carried 
out for a loop with x = 0.5 (filled down-triangles in Fig. ^). For low L, A$ w 1.8 $o? then 
A$(L) decreases substantially below the value for the filled disk, before increasing again, 
until the same period A$(L) is reached for L > 12. The loop is in the giant vortex state 
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when $ > 17 $o- For loops made of even thinner wires (x = 0.7 and x = 0.9) (open symbols 
in Fig. ||), T c ($) stays periodic up to L = 13. The constant period A$ corresponds to 
A$ m = ((1 + x)/2) 2 A$ = 1 $0; as it should be in the London limit for x — > 1, according 
to Eq. (0). 

In summary, we have solved the linearized GL equation for loops of different wire width, 
with Neumann boundary conditions at both the outer and the inner loop radius. The 
critical fields H* 3 (T) are always above the H C ^(T) = 1.69 H C 2(T), the surface critical field 
for a semi-infinite superconducting slab in contact with vacuum. The ratio H* 3 (T)/H c2 (T) 
increases when the size of the middle opening grows, i.e. in a sample topology with a large 
surface-to-volume ratio the nucleation field is strongly enhanced. Depending on the ratio 
inner to outer radius rj/r G of the loops, and on the applied magnetic flux, T c ($) shows a 
different behaviour: in thin-wire loops, the background of T c is parabolic (characteristic for 
2D behaviour) and the Little-Parks T c ($) oscillations are perfectly periodic. This behaviour 
can be described in the London limit. For loops with only a very small opening, the period 
of the T c ($) oscillations is decreasing with $ and the background T c reduction is quasi- 
linear (3D regime, e.g., like for a disk). Intermediate loops (for instance ri/r = 0.5) show 
a 2D-3D cross-over between the two regimes at a certain applied flux $ (corresponding to 
w = r —Ti = (1—x) r a w 2 £(T)), similar to the dimensional transition in thin films subjected 
to a parallel field. As soon the 3D regime is reached, a giant vortex state is created, where 
only a sheath close to the sample's outer interface is superconducting. The opening in the 
middle of the loop does not play a role anymore: T c ($) for the loop and for the disk become 
identical. 
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FIGURES 



FIG. 1. Calculated energy level scheme (dashed lines) for a superconducting loop with different 
ratio of inner to outer radius x = r.j/r : a) x = 0.1, b) x = 0.3, c) x = 0.5, d) x = 0.7, and 
e) x = 0.9. The lowest level for each magnetic flux $/$o corresponds to the highest possible 
temperature T C (H) for which superconductivity can exist. A state with phase winding number 
L = is formed at T c ($ ~ 0), and at each cusp in T C {H) the system makes a transition L — ► L + 1, 
indicating the entrance of an extra vortex. The solid and dotted lines correspond to H C 2(T) and 
H C 3(T), respectively. 

FIG. 2. Order parameter distribution 1^1 for x = ri/r Q = 0.1, at a fixed $ = 9<I>o- a) L = 2, 

b) L = 3, c) L = 4, d) L = 5, and e) L = 6. The latter (L = 6) corresponds to the ground state 
level | — ollll^ ^o) I; where the sample is in the giant vortex state. 

FIG. 3. Order parameter distribution |\&| for x = 0.5, at a fixed $ = 9 < I ) o- a) L = 2, b) L = 3, 

c) L = 4, d) L = 5, and e) L = 6. The state with L = 5 is only slightly modulated (l^j constant: 
Pa = 1-01) and corresponds to the ground state level | — ceLLL(9 &o)\- 

FIG. 4. Inverse enhancement factor rj" 1 = e(H* 3 ) / e(H C 2) for loops with different aspect ratio, 
compared to the case of a disk and an antidot. The horizontal dashed line at r/ _1 = 0.59 = 1/1.69 
corresponds to H c % (T) /H C 2 (T) = 1.69 for a plane superconductor /vacuum boundary. The insets 
show the supercurrent J profiles for the ground states at $ = 9$o f° r x = r i/ r o = 0.1 (left) and 
x = 0.5 (right). These are obtained from Eq. (^) and are normalized to —1 at the outer radius 
r = r Q . The corresponding |^| profiles are plotted in Fig. |2|e, and Fig. ||d, respectively. 

FIG. 5. The period A<& of the phase boundary T C (<I>) in units of the flux quantum <l?o as a 
function of the phase winding number L. The data, for several of the loops is shown as a symbol, 
and is compared to the period in a filled disk (filled square). The interconnecting lines are only a 
guide to the eye. 
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